REEXAMINATION OF SEVEN-DIMENSIONAL SITE PERCOLATION THRESHOLDS 



Dietrich Stauffer^ and Robert M. Ziff^ 

^ Institute for Theoretical Physics, Cologne University, D-50923 Koln, Euroland 
^ Department of Chemical Engineering, University of Michigan, Ann Arbor, MI 48109- 
2136, USA 

e-mail: stauffer@thp.uni-koeln.de, rziff@umich.edu 

Monte Carlo simulations alone could not clarify the corrections to scaling for the size- 
dependent Pc{L) above the upper critical dimension. Including the previous series estimate 
for the bulk threshold ]5c(oo) gives preference for the complicated corrections predicted by 
renormalization group and against the simple 1/L extrapolation. Additional Monte-Carlo 
simulations using the Leath method corroborate the series result for pc- 
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How to extrapolate from finite samples of linear dimension L to the thermodynamic 
limit L — i> cxD is an important problem in physics simulations. Above the upper critical 
dimension (4 in usual Ising models, 6 in random percolation), the critical exponents are 
known but nevertheless controversies remain^ for the five-dimensional Ising model. Earlier, 
for site percolation in seven dimensions, the numerical variation of the apparent threshold 
Pc{,L) = (pc) on a hypercubic lattice with sites could be fitted^ on Pc{oo) —pc{L) oc 1/L, 
but also'^ on the theoretical prediction Pc{oo) —pc{L) oc 1/L^ — const/L^/"^. Also for three- 
dimensional self-avoiding walks one may assume one empirical correction term, or many 
terms as predicted by renormalization methods with more free parameters^. Therefore 
the present note reexamines seven-dimensional site percolation in the search of a clarify- 
ing example, using a much better Cray-T3E computer and a slightly improved program^ 
compared with ref.2. 

We checked if the top of the hypercube is connected with the bottom hyperplane, using 
helical boundary conditions in five and free boundaries in the two remaining directions^ and 
a Hoshen-Kopelman algorithm without recycling^. Random number generation by integer 
multiplication with 16807 gave problems for L = 4, 8, 16 even for 64-bit arithmetic, but 
otherwise agreed with results from the Kirkpatrick-Stoll R250 generator; only the latter 
method was used for the data shown here. About 640 samples were averaged over for each 
L, giving an accuracy of about 10~^ for pc in large lattices. {[< pi > — < Pc >^]^/^ was 
not investigated in detail since it seems to vary as l/L"^.) Fig.l shows our results for L up 
to 20, together with a linear fit 

Pc{L) = 0.0909 - 0.12/ L (1) 

and the same fit as published before^: 

= 0.0887 - 3.3/L2 + 5/L^/2 . (2) 

We see that for large L both fits agree nicely with our new Monte Carlo results: On the 
basis of these simulations alone no preference for one fit over the other is visible. However, 
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when we extract from series expansions^ (in terms of reciprocal dimensionality) a bulk 
of about 0.089, then the more complicated fit (2) is clearly preferred over the linear fit (1). 
This is a nice example how the combination of series with Monte Carlo techniques can 
give some answers which are not clear from each technique separately^. However, a more 
direct and new series expansion determination of pc would be desirable. 

If we would have simulated L = 40 instead of L = 20 we might have chosen between 
the two fits on the basis of simulations alone. Such a size would have required about two 
orders of magnitude more in computer time and memory than the about 4000 processor 
hours (512 Mbyte each) used here. 

(About boundary conditions: If we check whether top and bottom are connected we 
cannot use vertical periodic boundary conditions. As long as not all directions are treated 
with periodic boundary conditions, and at least one direction uses free boundaries, the 
finite-size effects are expected to be dominated by the free boundaries. The bulk pc value 
is independent of the boundary conditions, and so may be the powers of L in finite-size 
corrections, but finite-size amplitudes as well as the fraction of samples spanning at Pc 
depend on such details.^ The Hoshen-Kopelman algorithm stores only one line of a square 
lattice at any one moment, and then it is simplest to use free boundaries also in this 
direction. In higher dimensions it is practical to store the L'^~^ sites of the one hyperplane 
kept in memory by a one-dimensional index i — 1, 2, ... , L'^^^ . In this way ^ we come to 
our mixture of helical and free boundary conditions.) 

Thus, to show that Monte Carlo can compete with series expansions, we made a 
different determination of the bulk pc using the Leath cluster growth algorithm^'^. Here, 
when a cluster stops growing before touching the boundaries, its properties are completely 
free of finite-size effects, and such data cannot be directly compared with the above Hoshen- 
Kopelman data for Pc{L). However, the bulk pc must be the same. We determined Pc by 
the methods given recently for three-dimensional lattices^^. We considered a virtual lattice 
of size 32^ and a maximum size cut-off of 16384 sites, and simulated 10^ - 10^ clusters 
each at various values of p. The clusters are grown independently for each value of p, since 
in this method the random numbers are not assigned to all lattice sites but only to those 
the cluster visits. 

It turned out that some of the largest clusters clusters (of maximum size of 16384 sites) 
wrapped a bit around the periodic boundary of the system, which could conceivably lead 
to the cluster touching itself and cause a bias. However, because of the very large number 
of sites in the lattice, the probability of this occurring is very low. (We could not easily 
check for wraparound error in our program.) In any case, the data for smaller clusters 
(where no wraparound was even possible) and for s — 16485 was completely consistent. 

Near the critical point, P>s, defined as the probability that a cluster grows to a size 
greater than or equal to s, behaves as 

P>s ~ s^- V((p - Pc)s") ~ s'-^iA + B{p - p,)s'' + . . .) (3) 

where r = 5/2 and a = 1/2 for d > 6. Thus in Fig. 2 we plot s^/^P>s vs. s^/^ for 
p = 0.0885, 0.0888, 0.0889 and 0.0890, and find good agreement with the expected behavior 
with A Ri 1.44, B ^ 20, and 

Pc = 0.08893 ± 0.00002 (4) 
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which is in exceUent agreement with the series extrapolation^. With this value, the simple 
fit of eq (1) can be excluded since it requires Pc = 0.091. 
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Fig.l: Effective thresholds PcXL) versus 1/L, together with eqs.(l,2). Part a gives the 
overall trend, part b expands the region of large lattices. The horizontal line is the series 
extrapolation for infinite L. 

Fig.2: Plots of s^/^P>s vs s ^j^^ Lg^th simulations, for p = 0.0890, 0.0889, 0.0888, 
and 0.0885 from top to bottom. The slopes of the three upper curves are 0.0014, -0.0006 
and -0.0026 respectively. 



3 



0.095 



0.09 - 



0.085 - 



0.08 - 



0.075 - 



0.07 - 



0.065 







Effective site percolation tliresliolds for L'^7 site, and series result 

T 




0.05 



0.1 



0.15 



0.2 



0.25 



1/L 



T 



Expanded central region of this figure 
1 1 



0.085 - 



+ 



0.084 - 



o 



0.083 - 



0.082 - 



0.081 - 



0.08 - 



0.05 



0.06 



0.07 



1/L 



0.08 



0.09 



0.1 



